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Abstract 

We discuss single heavy neutrino production e+e“ —> Nn —> iWn, i = e, g, r, 
at a future high energy collider like CLIC, with a centre of mass energy of 3 TeV. 
This process could allow to detect heavy neutrinos with masses of 1—2 TeV if their 
coupling to the electron VeN is in the range 0.004—0.01. We study the dependence 
of the limits on the heavy neutrino mass and emphasise the crucial role of lepton 
flavour in the discovery of a positive signal at CLIC energy. We present strategies 
to determine heavy neutrino properties once they are discovered, namely their 
Dirac or Majorana character and the size and chirality of their charged current 
couplings. Conversely, if no signal is found, the bound VeN < 0.002 — 0.006 would 
be set for masses of 1 — 2 TeV, improving the present limit up to a factor of 30. We 
also extend previous work examining in detail the flavour and mass dependence 
of the corresponding limits at ILC, as well as the determination of heavy neutrino 
properties if they are discovered at this collider. 


1 Introduction 

The existence of heavy neutrinos is usually associated to the see-saw mechanism [1], 
which provides a simple and elegant explanation for the smallness of light neutrino 
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masses. This economical solution has no phenomenological implications at large collid¬ 
ers, however, because the new neutrinos are extremely heavy, with masses of the order 
of 10^^ GeV for Yukawa couplings Y of order one. These extra neutrinos also sup¬ 
ply a mechanism to explain the observed baryon asymmetry of the universe through 
leptogenesis [2]. Many attempts have been made to construct viable models of this 
type but with neutrino masses at the TeV scale [3,4]. The price to pay in all cases 
is the loss of simplicity. Heavy neutrinos give contributions to light neutrino masses 
of the order Y'^v'^/mN, where v is the vacuum expectation value of the Higgs boson 
and niN the heavy neutrino mass. These contributions are far too large for iutv in the 
TeV range unless (i) Y is very small, of order 10“^, in which case the heavy neutrino 
is almost decoupled from the rest of the fermions, or (ii) there is another source for 
neutrino masses giving a comparable contribution cancelling the ~ Y‘^v ‘^/ttin one from 
the see-saw mechanism. 

Both solutions require a theoretical effort so as to build a natural model which 
reproduces light neutrino masses. In the hrst case, it is necessary to justify why neutrino 
Yukawa couplings are much smaller than for charged leptons and quarks. In the latter, 
not only it is necessary to provide an additional source of neutrino masses, but it is also 
crucial to give a natural explanation for the (apparently hne-tuned) cancellation of both 
contributions [5]. But despite the disadvantage of complexity there is the signihcant 
beneht that these models, experimentally not excluded, might be directly testable at 
future colliders by searching for the production of heavy neutrinos. (Additionally, there 
could be indirect evidence of their presence in neutrino oscillation experiments [6].) An 
important question is then whether these heavy states are indeed observable or not. 
Although their masses are within the reach of forthcoming or planned colliders, their 
mixing with the Standard Model (SM) leptons must be also large enough to allow for 
their production at detectable rates. This is because they are SM singlets, and in the 
absence of new interactions their couplings are proportional to their mixing with the 
light neutrinos. 

Independently of the mass generation mechanism, heavy neutrinos with masses of 
several hundreds of GeV appear in Grand Unihed Theories, like for instance in those 
based on SO(IO) or on larger groups as E(6) [7], and can survive to low energies 
[8]. Kaluza-Klein towers of neutrinos are also predicted in models with large extra 
dimensions, being possible to have the hrst heavy modes near the electroweak scale [9]. 
Their existence is allowed by low energy data, which set strong constraints on their 
mixing with the light leptons but leave room for their production and discovery at 
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large colliders. If they have masses up to 400 GeV and a mixing with the electron 
VeN ~ 0.01, they will be discovered at an international linear collider (ILC) with a 
centre of mass (CM) energy of 500 GeV [10]. An eventual ILC upgrade to 800 GeV will 
extend the reach to higher masses, but in order to experimentally test the existence 
of TeV scale neutrinos a larger CM energy is required, which is achievable only at a 
future collider in the multi-TeV range, like the compact linear collider (CLIC) 

with a CM energy of 3 TeV [11,12]. 

In this paper we present a study of the CLIC potential to discover heavy neutrino 
singlets and determine their properties in the process > Ni/ —> IWv. In section 

[2] we review the formalism and derive the interactions of heavy Dirac and Majorana 
neutrino singlets with the gauge and Higgs bosons, summarising present constraints on 
their couplings to the charged leptons. In sectionISjwe discuss the general characteristics 
of elTz/, jjWv and rWu final states. We analyse the different contributions to the signal 
and background, stressing the crucial fact that an eNW coupling is necessary to observe 
the heavy neutrino in any of the channels. In section H) we describe the procedure 
used for our Monte Carlo calculations. The sensitivity and limits on charged current 
couplings achieved at CLIC are discussed in section 0 examining also the dependence 
on m^v. In case that a heavy neutrino was discovered, we show how its Dirac or 
Majorana nature could be established and its charged couplings measured. In section 
El we perform a similar analysis for neutrinos with masses of 200 — 400 GeV at ILC, 
extending previous work [10], also studying the determination of their properties and 
comparing with CLIC results. In section Owe draw our conclusions. 

2 Addition of neutrino singlets 

In this paper we consider a SM extension with heavy Majorana (M) or Dirac (D) 
neutrino singlets. The most common situation is that three additional heavy eigenstates 
Ni, i = 1, 2, 3 are introduced, and for definiteness this is what we assume in this section. 
The formalism is however general for any number of singlets [14]. In the following we 
will obtain the interactions of heavy neutrinos with the light leptons, pointing out the 
differences between the Dirac and Majorana cases when they exist. 

The neutrino weak isospin T3 = 1/2 eigenstates ^'iL are the same as in the SM. In 
the case of Dirac neutrinos we introduce 9 additional SU(2)2, singlet fields 

Al ) ^iRi^iRi * = 1)2,3, (D) (1) 
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which allow the light neutrinos to have Dirac masses too. For Majorana neutrinos only 
three helds are added 

N'^, t = 1,2,3, (M) (2) 

with = iyliY, NY = {NYY- In matrix notation, the form of the mass terms in the 
Lagrangian is similar in both cases, 

T) (I") 

^ i;;) (;«)+H.C,. (M) ,3) 

where the Y, B and M blocks are 3x3 matrices.^ The physical meaning is of course 
different, since the Y matrices correspond to Yukawa interactions, B, B' are bare mass 
terms and Ml, Mr lepton number violating Majorana mass matrices.^ The complete 
mass matrices M. can be diagonalised by U^j^M.Ur = A^diag- For Majorana neutrinos 
IAr = Ur, while for Dirac neutrinos the two unitary matrices are independent. The 
mass eigenstates are 



Both for Majorana and Dirac neutrinos the weak interaction Lagrangian is written in 
the weak eigenstate basis as 

Cw = -^47XIF^ + H.c., 

= -^v'lI^v'lZ^,, (5) 

^Both mass terms in Eqs. m are particular cases of a general 12 x 12 symmetric mass matrix 
connecting (Y YrY N'^ and {{vrY {^lY Yjj). If we assign lepton numbers L = 1 to 

L = —1 to ,{NYY, the mass term for Dirac neutrinos corresponds to the L = 0 
entries (with conserved lepton number). The mass term for Majorana neutrinos includes the Yukawa 
entries Y with L = 0 and the diagonal lepton number violating blocks Mr, Mr with L = 2, L = —2, 
respectively, and in this case i/jj, N'j^ are assumed very heavy or decoupled. 

^In the Dirac case the right-handed states are equivalent, and by a suitable redefinition 

one can always choose a weak basis with B' = 0. Additionally, with adequate rotations B and Mr 
could be assumed diagonal without loss of generality. This is not necessary for our discussion, anyway, 
and the results in this section do not rely on such assumptions. 
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with I'^j^ the charged lepton weak eigenstates. Let ns divide for convenience the rotation 
matrices Ul, Ur in 3 x 6 blocks, 


Ul = 


Ul 

K 




( 6 ) 


Then, the weak interaction Lagrangian is written in the mass eigenstate basis as 


Cw 





2o 


w 


[Vl Nl) Y uIUl 


+ H.c., 



(7) 

( 8 ) 


where t/; is a 3 x 3 nnitary matrix resnlting from the diagonalisation of the charged 
lepton mass matrix. The extended Maki-Nakagawa-Sakata (MNS) matrix [13] V = 
uJUl has dimension 3x6. Neutral interactions are described by the 6x6 matrix 
X = U\Ul, related to the former by X = V^V. 

The interactions with the Higgs boson II are 


Hh — — 


I^H — 


1 

7! 

1 

7! 

1 

V2 


{v'^YN'^ + H + H.c., 

{VlNl)uI{YU'r + Y'Ur) ' 


X 


H + H.c., (D) 


R 


Yj^YN^HY H.C. 


^R 




(M) 


(9) 


In order to obtain explicit expressions in terms of masses and mixing angles, we decom¬ 
pose V in two 3x3 blocks, V = with parameterising the mixing 

of the charged leptons with the light and heavy neutrinos, respectively. The latter is 
experimentally constrained to be small (see below), thus terms of order can 

be neglected. After a little algebra, the scalar interactions of both heavy Dirac and 
Majorana neutrinos can be written as 

= -YYT + H.c., (10) 

Z1V1\Y 

with Mjv their 3x3 diagonal mass matrix Mjv = diag (mAr^, ttiatj, mAfg). In the Dirac 
case there are additional Yukawa couplings among the light neutrinos. 
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The mixing of heavy neutrinos with charged leptons is restricted by two groups of 
processes [15-21]: (i) tt ^ £iy, Z ^ ui/ and other tree-level processes involving light 
neutrinos in the final state; (h) /r —*• ey, Z —> £~^£'~ and other lepton flavour violating 
(LFV) processes to which heavy neutrinos can contribute at one loop level. All these 
processes constrain the quantities 

3 3 

^ yivyi'ui = X] ^eNiVgl]^.. (11) 

i=l i=l 

The processes in the first group measure lepton charged current couplings. A global 
fit yields the bounds [20] 

Qee < 0.0054 , < 0.0096 , Qrr < 0.016 , (12) 

with a 90% confidence level (CL). For heavy neutrino masses in the TeV range, LFV 
processes in the second group give the constraints [19] 

< 0.0001, |fler|<0.01, |fl^^|<0.01. (13) 

The limits in Eqs. ca are model-independent to a large extent, and independent of 
heavy neutrino masses as well. They imply that the mixing of the heavy eigenstates 
with the charged leptons is very small, — 0.0054, 0.0096, 0.016 for £ = e, fi, r, 

respectively. On the other hand, the bounds in Eqs. ina do not directly constrain the 
products ViNiV^iNi Out the sums in the r.h.s. of Eq. (fTTll . and cancellations might 
occur between two or more terms, and also with other new physics contributions. 
These cancellations may be more or less natural, but in any case such possibility 
makes the limits in Eqs. 03 relatively weak if more than one heavy neutrino exists 
[10, 22]. Besides, we note that these limits are independent of the heavy neutrino 
nature. Eor heavy Majorana neutrinos there is an additional restriction from the non¬ 
observation of neutrinoless double beta decay, which is below present experimental 
limits for |14ivP < 0.0054 and > 100 GeV [23]. 

Since mixing of the charged leptons with heavy neutrinos is experimentally required 
to be very small, the usual MNS matrix is approximately unitary, up to corrections 
of order Moreover, at large collider energies the light neutrino masses can be 

neglected. With these approximations can be taken equal to the identity matrix, 
implying also = 6u', = Wwi, be. the vertices between light leptons can be 

taken equal to their SM values for massless neutrinos, and the couplings for flavour- 
changing neutral interactions v^NiZ are proportional to those for charged currents 
£NiW. 
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The production of a heavy neutrino N involves its interactions with the light 
fermions. The charged current vertex with a charged lepton i can be directly read 
from Eq. o. 

Cw = {iYVmPLN + NYV;„PLe Wl) . ( 14 ) 

The neutral current gauge couplings with a light neutrino are 

Cz = {ua^VeNPiN + nyv;^Pli^,) . (15) 

ZC\Y 

In the Dirac case, the two terms in Cz describe the interactions of heavy neutrinos 
and antineutrinos. If they are Majorana particles, the second term can be rewritten in 
terms of and N, giving 

Cz = -7^m^iViNPL-v;^Pn)N Z^. (M) (16) 

Zcw 

The scalar interactions of the heavy neutrino are 

VmPnN + N H , (17) 

ZMw 

where the second term can again be rewritten for Majorana neutrinos, 

= + H. (M) (18) 

ZM]y 


For our computations it is also necessary to know the total heavy neutrino width 
T^. N can decay in the channels N —> (if is a Majorana fermion N —> 

is allowed as well), N ^ Zuf^ and N —> Hui. The partial widths for these decays 
are [22,24] 


r(iv ^ lE+r) 

ToiN^Zue) 
> Zui) 
Tn{N ^ Hu,) 

TMiN ^ Hu,) 


r(iv 

9^ 


w-e^ 


Gdvr 


1 14. 


m 


'TV 


N 






1287rc^ 


1 14. 


'N\ 


12 ^N 

M| 


2 r,,(iv- 

; 114 AT 
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■ Zu ,), 

2^N 


1287r' 

2r^(iV 




Hu,), 


m 


N 




2^ 


1 - 


M| 


m 


N 


1 - 




m 


N 


m 


N 


1 + 


Ml 


m 


m 


N 


- 2 - 


Mi 


N 


m‘ 


N 


( 19 ) 
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with r£) and I'm standing for the widths of a Dirac and Majorana nentrino. The factors 
of two in the partial widths oi N ^ N ^ Hp^ for a Majorana nentrino are the 
consequence of the extra couplings in Eqs. Hi,!!!!), which are not present for a 
Dirac neutrino. From Eqs. (HID it follows that for equal values of the mixing angles 
ViN the width of a heavy Majorana neutrino is twice as large as for a Dirac neutrino. 
Another straightforward consequence of these expressions is that the partial widths for 
IT, Z and Higgs decays are in the ratios 2:1:1 (the latter for mn -C ruTv). Since the 
Higgs mass is still unknown, we will ignore the decays N —> Hp(^ in the calculation of 
Etv. If these decays are included, the branching ratios (and hence the hnal signal 

cross sections) are multiplied by a factor which ranges between 3/4 (for niH -C mAr) 
and unity (for niH > 


3 iWiz production and lepton flavour 

The existence of new heavy neutrinos is rather difficult to detect as an excess in the 
total cross section for iqq'p. The stringent experimental bounds on their 

mixing angles with the light particles restrict the size of their contribution to this 
process to a few percent, except for low rriN values. Such a small increase in the 
cross section is unobservable due to the inherent uncertainties in the SM prediction. 
Nevertheless, the heavy neutrino contribution to this signal is dominated by on-shell 
N production [22,25] e+e“ —> Np —> i~W^p —> £~qq'p (plus the charge conjugate) 
if kinematically accessible, yielding a peak in the iqq' invariant mass distribution. If 
heavy neutrino mass differences are of the order of 100 GeV or more the neutrino peaks 
do not overlap, so that their experimental study can be done independently, since in 
this case the interference of the relevant amplitudes is negligible. We will thus assume 
for simplicity that only one heavy neutrino N is produced. For quasidegenerate heavy 
neutrinos with -C 1 [4] one must consider interference effects, 

which are not addressed here. 

In Np production the two quarks in the hnal state result from the decay of an 
on-shell W boson, hence we can safely conhne the analysis to the phase space region 
where their invariant mass is not far from Mw, and restrict the calculation to £Wp 
production (with IT — qq') in the presence of a heavy neutrino. We hrst discuss the 
process for £ = e and later point out the differences for £ = p, r, using a reference value 
mTv = 1500 GeV. For most purposes £~ and production may be summed because 
CP-violating effects are negligible, as it is briehy commented at the end of this section. 
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3.1 Final states with electrons 


The Feynman diagrams for e’''e“ —> e~W^v involving heavy neutrino exchange are 
shown in Fig. [H We neglect the electron mass (as well as light quark masses) in the 
computation of matrix elements, thus we do not include scalar diagrams. When the 
outgoing light neutrino flavour is constrained to be Ve we write it explicitly. Diagrams 
[lb and m are present only if iV is a Majorana fermion, while the rest are common to 
the Dirac and Majorana cases. 




Figure 1: Feynman diagrams for e’''e —> e mediated by a heavy neutrino N. 

The differences between the Dirac and Majorana cases are discussed in the text. 

We notice that the six diagrams have an eNW vertex, whose presence is thus 
necessary for this process to occur. The three first ones involve the production of 
on-shell N, the corresponding amplitudes being proportional to Ww times a factor 
0(1) from the N decay branching ratio. Their contribution to the cross section is 
then quadratic in Ww (interference with the SM amplitude is negligible). The last 
three diagrams are proportional to Fjv, giving V^v terms in the cross section (through 
interference), plus ITjv terms. With present limits on Ww, the former are one order of 
magnitude larger than the latter, but still remain two orders of magnitude below the 
size of V^v terms from diagrams with N on its mass shell. Among these, the t- and 
M-channel W exchange diagrams ([lb and Hb, respectively) are the only ones relevant 
at CLIC energy and s-channel Z exchange m is highly suppressed, being a factor 
~ 5 smaller than off-shell contributions. We then arrive to one crucial point: for equal 
values of the mixing angles, the contributions of Majorana and Dirac heavy neutrinos 
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to the e“ggV cross section are almost equal. The reason is simple: for a Majorana N 
the neutrino signal is strongly dominated by two non-interfering Feynman diagrams 
which give equal contributions to the cross section, instead of one in the Dirac case. 
On the other hand, the width of a Majorana N is twice as large as for a Dirac one, as 
noted in the previous section. 


In the phase space region of interest, the relevant SM diagrams are those for 
e+e“ — e~W~^h'e with subsequent hadronic W decay, as depicted in Fig. El The 
main contribution comes from diagram El and is one order of magnitude larger than 
the rest. This fact constrasts with the behaviour at ILC energy, where about one half 
of the cross sections comes from resonant W~^W~ production, especially from diagram 
Ell. The 8 additional diagrams for e’''e —e qq'u which do not involve the decay 

W~^ qq' give an irrelevant contribution (smaller than 0.5 %) in the phase space 
region studied. The quadratic corrections to the iuW and vertices neglected in 
section El are unobservable with the available statistics. We also ignore diagrams like 
[Bd and m with the exchange of a light Majorana neutrino, which are suppressed by the 
ratio rriu/Q ~ 10“^^, where Q ~ ^/s is the typical energy scale in the process. 


It is worth pointing out the effect of beam polarisation on the cross sections for 
e~W^v production through N exchange only and through SM diagrams. For N ex¬ 


change we have a+- : cr„+„- = 44000 : 1 (for rriM = 


a^+- = a+- = 0. For the SM process, we find : 

^R^R ^R^L 


1500 GeV, Kjv = 0.05), with 
^e+e- : = 3^00 : 200 : 

^R^R ^L^R 


1, o'e+e- = 0. Therefore, the use of negative electron polarisation P^- and positive 
positron polarisation Pe+ improves the observability of the signal as well as the statis¬ 
tics. Besides, it can be seen that e^e)j cross sections only receive contributions from 
diagrams with Z or photon s-channel exchange, thus the huge difference between 
and ejejj cross sections reflects the suppression of the former. 


3.2 Final states with muons and tans 

Final states with £ = /r, r share similar production properties and we refer to muons for 
brevity. The diagrams for e’''e“ — jj,~W^iy via heavy neutrino exchange are shown in 
Fig. El The same comments regarding the contributions of the diagrams in Fig. [B apply 
in this case (up to different mixing angles). We observe that all contributions except Et 
involve an electron-heavy neutrino interaction; in particular, the leading diagrams Eli 
andEb correspond to heavy neutrino production via an eNW vertex with subsequent 
decay through a /j,NW interaction. This leads to the important consequence that the 
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Figure 2: Feynman diagrams for e'''e —^ e within the SM. 


signal of heavy nentrinos is relevant only if N simnltaneonsly mixes with the 
electron and mnon. We also notice that, even withont mixing with the mnon, a heavy 
neutrino can mediate production, via diagram Etl. Nevertheless, the cross 

section is very small in this case. The SM background is with 6 

Feynman diagrams like the ones in[2K-|2lF bnt replacing e“, z/g by /i“, z/^, respectively. 
Its cross section is dominated by resonant W~^W~ prodnction and is 30 times smaller 
than for e~W^Ue. 

In order to discuss the effect of beam polarisation we have to distingnish two cases. If 
prodnction takes place mainly throngh a eNW coupling, beam polarisation has 
clearly the same effect as in e~W~^h' production: negative P^- and positive Pe+ enhance 
the signal and background and thus improve the statistics. If the heavy neutrino does 
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w+ 
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Figure 3: Feynman diagrams for e'''e —fi mediated by a heavy nentrino N. 

The differences between the Dirac and Majorana cases are discussed in the text. 

not mix with the electron bnt mixes with the mnon, the reverse sitnation occnrs. 
Since the only contribntion comes from diagram Et, the nse of left-handed positrons 
and right-handed electrons actually increases the signal, while reducing the SM cross 
section for this process [10]. This case has no interest at CLIC energy, anyway, because 
for VeN = 0 the signal is not observable. 

3.3 CP violation with heavy neutrinos 

In this work we do not address any CP violation effects, which are unobservable in 
the processes studied. Partial rate asymmetries between £~ and final states are 
negligible at the tree level. They reqnire interference of diagrams with different CP- 
conserving phases {e.g. diagram^, with a phase arising from the N propagator [27], 
and one of the diagrams in Fig. H- This interference is very small due to kinematics 
(the width Fat is small compared to the energy scale). Another possibility is the study 
of triple-product asymmetries in the decay of the heavy nentrino. However, these 
asymmetries are proportional to the mass of the final state charged lepton [28], and 
hence very small. Therefore, we snm i~ and final states in all onr results, unless 
otherwise stated. 
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4 Event generation 


The matrix elements for > i~qq'u are calculated using HELAS [29], 

including all spin correlations and finite width effects. We sum SM and heavy neutrino- 
mediated diagrams at the amplitude level. For Majorana fermions the Feynman rules 
are given in Ref. [30]. We assume a CM energy of 3 TeV, with a beam spread of 0.35% 
[11], ignoring an eventual beam crossing angle of 0.02 rad at the interaction point [12] 
whose effect in our simulation is very small. We make use of electron polarisation 
Pe- = —0.8 and positron polarisation Pe+ = 0.6. The luminosity is taken as 1000 fb“^ 
per year. In order to take into account the effect of initial state radiation (ISR) and 
beamstrahlung we convolute the differential cross section with “structure functions” 
Disnix), D^six), 


da = da{xiyiE, 0 : 21 / 2 ^) AsR(a^i)^Bs(f/i)AsR(a^ 2 )^Bs(f/ 2 ) dxidyidx 2 dy 2 ■ 

Jo 

The function describing the effect of ISR used is [31] 

Disk{x) = ^( 1 - a ;)^ 


( 20 ) 


4(1-7) 


X 


r(i + f) 

1(1 + - I (1(1 + 3i^) logi - (1 - xY 


where 


p(s) = —6 log 


, Oio s 

1 - ^ log — 

37 r mi 


( 21 ) 


( 22 ) 


7 is the Euler constant, ao = 1/137 the fine structure constant, s the center of mass 
energy squared and me the electron mass. For beamstrahlung we use [32] 


^Bs(a:) = e 


-N 


„—k{1—x)Ix 


/I \ 1 , (23) 

a;(l — X) J 

with N = N^/2, K = 2/3 T. For the proposed luminosity we take the parameters 
T = 8.1, = 2.3 [11]. The function h{y) is [33] 




^ n\ F(n/3) ’ 

where y = N[k (1 — x)/xY^^. For large y, h{y) has the asymptotic expansion 

35 1295 


(24) 




X 16588x2 


+ 


(25) 
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with 2; = (|//3)^/"^. 

In final states with r leptons, we select r decays to n, p and ai mesons (with a 
combined branching fraction of 54% [34]), in which a single Ur is produced, discarding 
other hadronic and leptonic decays. We simulate the r decay assuming that the meson 
and r momenta are collinear (what is a good approximation for high r energies) and 
assigning a random fraction x of the r momentum to the meson, according to the 
left-handed probability distributions [35] 

P{x) = 2{l-x) (26) 


for pious, and 

= 2(3 -4(2 + 1 [(1 - 2C") - (1 - 203^] (27) 

for p and Oi mesons, where C = ^p,ai/^T- We assume a r jet tagging efficiency of 
50%. In certain measurements the use of c tagging is necessary as well, and we assume 
a 50% efficiency, the same one that it is expected at ILC [36]. 

We simulate the calorimeter and tracking resolution of the detector by performing 
a Gaussian smearing of the energies of electrons (e), muons (/i) and jets (j), using the 
expected resolutions [12], 




10 % 


© 1 % , 


AE^^ 

Ef^ 


0.005% , 


AE^ 



(28) 


where the two terms are added in quadrature and the energies are in GeV. We apply 
kinematical cuts on transverse momenta, pr > 10 GeV, and pseudorapidities |?7| < 2.5, 
the latter corresponding to polar angles 10° < 9 < 170°. We reject events in which the 
leptons or jets are not isolated, requiring a “lego-plot” separation AR = > 

0.4. For the Monte Garlo integration in 6-body phase space we use RAMBO [37]. 

In electron and muon final states the light neutrino momentum pj, is determined 
from the missing transverse and longitudinal momentum of the event and the require¬ 
ment that pI = 0 (despite ISR and beamstrahlung, the missing longitudinal momentum 
approximates with a reasonable accuracy the original neutrino momentum). In final 
states with r leptons, the reconstruction is more involved, due to the secondary neu¬ 
trino from the r decay. We determine the “primary” neutrino momentum and the 
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fraction x of the r momentum retained by the r jet using the kinematical constraints 

E]y + H —Ej = \/s , 

X 

^ 1 ^ 

Pw+Pv + -Pj = 0 , 

X 

pI = 0, (29) 

in obvious notation. These constraints only hold if ISR and beamstrahlung are ignored, 
and in the limit of perfect detector resolution. When solving them for the generated 
Monte Carlo events we sometimes obtain a; > 1 or a: < 0. In the hrst case we arbitrarily 
set a; = 1, and in the second case we set x = 0.55, which is the average momentum 
fraction of the r jets. With the procedure outlined here, the reconstructed r momen¬ 
tum reproduces with a fair accuracy the original one, while the obtained is often 
completely different from its actual value. 


5 Heavy neutrino discovery at CLIC and determina¬ 
tion of their properties 

We address in turn the discovery of a new heavy neutrino fsections 15.11 and 15.211 . the 
determination of its Dirac or Majorana character (section 15.511 and the measurement 
of its couplings to e, p, r fsection . We try to be as concise as possible without 

losing generality or omitting the main points. Following the discussion in section |H1 
we can distinguish two interesting scenarios for our analysis: (i) the heavy neutrino 
only mixes with the electron; (n) it mixes with e and either p, r, or both. For the 
study of the niN dependence and the determination of the neutrino nature we assume 
for simplicity that N only mixes with the electron. Additionally, we assume that the 
neutrino is a Majorana fermion in sections lOl 15.21 and 15.41 where the results obtained 
are almost independent of its character. 

5.1 Discovery of a heavy neutrino 

The existence of a heavy neutrino which couples to the electron can be detected as 
a sharp peak in the distribution of the ejj invariant mass mejj, plotted in Fig. 21 for 
VeN = 0.05. The dotted and solid lines correspond to the SM and SM plus a 1500 GeV 
Majorana neutrino, respectively. For a Dirac neutrino the results do not change. The 
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width of the peak is mainly due to energy smearing included in our Monte Carlo, and 
the intrinsic N width, Fat = 8.2 GeV and Fjv = 4.1 GeV for a Majorana and Dirac 
neutrino, respectively, has a smaller influence in this case. 



Figure 4: Kinematical distribution of the ejj invariant mass, for ruAr = 1500 GeV. 


The SM and SM plus heavy neutrino cross sections are collected in Table [H before 
and after the kinematical cut 

1460 GeV < rriejj < 1540 GeV. (30) 

The criterion used here for the discovery of the new neutrino is that the excess of 
events^ (the signal S) in the peak region dehned by Eq. pUj) amounts to more than 
5 standard deviations of the number of expected events (the background B), that is, 
S/y/B > 5. This ratio is larger than 5 for V^n > 7.8 x 10“^, which is the minimum 
mixing angle for which a 1500 GeV neutrino can be discovered. Conversely, if no signal 
is found, the limit VeN < 4.5 x 10“^ can be set at 90% conhdence level (CL), improving 
the present limit VeN < 0.073 by a factor of 16. 

No cut With cut 
SM 516 14.6 

SM + N 548 39.4 

Table 1: Cross sections (in fb) for before and after the kinematical 

cut in Eq. (EDI). 


^It must be stressed that the SM cross section at the peak can be calculated and normalised using 
the measurements far from this region. 


16 












In the most general case that N simultaneonsly mixes with the three charged lep¬ 
tons, there may be in principle signals in the e, fi and r channels, and the three of them 
must be experimentally analysed. We choose equal values Ww = Vfj,N = Wiv = 0.04 
to illustrate the relative sensitivities of the three channels. For electron and muon 
hnal states we apply the kinematical cut in Eq. (HD, while for tans the distribution is 
broader and we use 

1420 GeV < nirjj < 1580 GeV. (31) 

The SM and SM plus heavy neutrino cross sections after these kinematical cuts"^ can 
be found in Table El For these values of the couplings, the heavy neutrino signal could 
be seen with a statistical signihcance of ~ 40a, ~ 250a and ~ 70cr in the e, fi and r 
channels, respectively, after one year of running. 

e fi T 

SM 14.6 0.36 0.096 

SM + N 19.5 5.24 1.19 

Table 2: Cross sections (in fb) for for i = e,fi,T, including the 

kinematical cuts on mijj. 


We clearly see that the muon channel is much more sensitive for equal values of the 
couplings. An eNW interaction is absolutely necessary to produce the heavy neutrino 
at observable rates but, once it has been produced, its decays in the muon and tan 
channels are easier to spot over the background. Since the production mechanism 
is strongly dominated by t and u-channel W exchange diagrams (see section H, the 
observed signals Se, S^j,, Sr can be written as 


S, = ApV 


^eN 


T/2 


v:-„+Vi 


liN 




N 


(32) 


to an excellent approximation. The common factor 1/^ comes from the production, 
the ratio of couplings corresponds to the decay and are constants. Using the data 
in Table O and Eqs. (El we can obtain the combined limits on VeN and or U-tv 
plotted in Fig. El It is very interesting to observe that a small coupling to the muon 

^For the muon and tau channels the background is dominated by resonant W~^W~ production, 
thus a cut on the invariant mass could reduce it significantly. However, in practice it may be very 
difficult to reconstruct the W mass at CLIC energy, and to be conservative we do not apply any cut 
on 
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V^j_N ^ 0.005 greatly increases the sensitivity to Kat, from ~ 0.008 to ~ 0.0035, dne 
to the better observation of the heavy nentrino in the mnon channel. For a fixed V(.n 
and increasing V^tv, the eWv channel becomes less significant because of the smaller 
branching ratio, and the observation is better in the v channel. (For ILC the 
behaviour is rather different, see next section.) In the case of the tan the effect is 
similar but less pronounced. 



0.01 


(a) 


(b) 


Figure 5: Combined limits on: V^n and for VrN = 0 (a); V^n and Ktv, for = 0 
(b), for rriAT = 1500 GeV. The red areas represent the 90% CL limits if no signal is 
observed. The white areas extend up to present bounds VeN < 0.073, I%7v < 0.098, 
Vtn < 0.13, and correspond to the region where a combined statistical significance of 
5a or larger is achieved. The indirect limit from p. — e LFV processes is also shown. 


For comparison, we include in Fig. El (a) the indirect limit on I4 jv, VfiN derived 
from low energy LFV processes. For V^n smaller than 0.05, the direct limit obtained 
by the absence of heavy neutrino production is much better than the indirect one. 
Moreover, the latter can be evaded if we allow for cancellations between heavy neutrino 
contributions, as discussed in section [2l (see also Ref. [10]). In the case of the tan 
lepton, the indirect limit from LFV processes is less stringent than the direct limits 
VeN < 0.073, VrN < 0.13 and is not shown. 

When the heavy neutrino does not couple to the electron but only to the muon 
and/or tan, the only Feynman diagram contributing to the signal isOt. For CLIC the 
situation is much worse than for ILC [10] because at a higher CM energy this diagram 
is more suppressed. Heavy neutrino production rates are thus negligibly small, giving 
an excess of a handful of events (for an integrated luminosity of 1000 fb~^) over the 
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SM background even for V^at, VrN in their upper experimental bounds. 

5.2 Dependence on the heavy neutrino mass 

The cross section for e+e“ —including the heavy neutrino contribution ex¬ 
hibits a moderate dependence on ruAr, which is mainly due to phase space. The vari¬ 
ation of the total cross section (including ISR, beamstrahlung and “detector 

cuts” as explained in section with mjq can be seen in Fig. ini(a). For larger masses 
the cross sections are smaller an thus the limits on VeN are worse. This behaviour is 
attenuated by the fact that the SM background also decreases for larger mejj, as can 
be clearly observed in Fig. HI The resulting limits on the heavy neutrino coupling are 
shown in Fig.|ni(b) as a function of mAr. These limits assume implicitly that the heavy 
neutrino only mixes with the electron (for mixing also with the muon they improve, 
as seen in the previous subsection). The kinematical cuts are not optimised for each 
value of mAT. 

700 


650 

e 

^ 600 

550 


1000 1500 2000 2500 3000 1000 1500 2000 2500 3000 

% 

(a) (b) 

Figure 6: (a) Cross section for e+e" —^ for V^n = 0.05 and different values of 

rriN- (b) Dependence of the discovery and upper limits on VeN on the heavy neutrino 
mass. Both plots assume mixing only with the electron. 




5.3 Determination of the Majorana or Dirac nature 

The cross section for e~W^v production mediated by a heavy neutrino is fairly insen¬ 
sitive to its Dirac or Majorana character. Still, the different production mechanisms 
show up in the angular distribution of N with respect to the incoming electron. In 
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the case of a heavy Majorana neutrino the new contribution is dominated by diagrams 
[T^, mb, leading to a forward-backward symmetric distribution for the production angle 
(pNe- between N and the incoming electron. The distribution peaks at cos ip^e- = 1 
when t = {pn — Pe-Y = 0 and the hrst diagram is enhanced, and at cosipNe- = —1 
when u = {pN — Pe+)‘^ = 0 (and the second one is enhanced). 


In the case of Dirac neutrinos the u-channel diagram is absent, and the distribution 
only peaks at cosip]\fe- = 1 for hnal states with e~,p~,T~ and at cosp^e- = 
for the charge conjugate processes with e~^, . It is then convenient to dehne 

= PNe-T‘^Ne+ for e~W^h', e'^W~h' events, respectively. Its normalised distribution 
is shown in Fig.m(a) for events surviving the kinematical cut in Eq. jHi. We consider 
the SM, and SM plus a Majorana or Dirac neutrino. The most conspicuous results are 
obtained subtracting the SM contribution, which can be calculated and calibrated using 
the measurements outside the peak. In this case. Fig. m(b), the difference between 
the Dirac and Majorana cases is apparent. We remark that the signal cross section at 
the peak is of 24.8 fb for VeN = 0.05, what gives a sufficiently large event sample to 
distinguish both cases even for smaller mixing angles (see also Fig. in])- 



COS(Ppj 


1 


0.8 



0.2 


— N(M) 
-- N(D) 


0 

cos (p 
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(a) 


(b) 


Figure 7: (a) Dependence of the cross section on the angle pn, for the SM and the SM 
plus a 1500 GeV Majorana (M) or Dirac (D) neutrino, (b) The same, but subtracting 
the SM contribution. 


For completeness, we also show in Fig. [HI the dependence of the cross section on 
the angles pw, Pe between the produced and the incoming electron (for e~ 

hnal states) or positron (for hnal e"*"). We restrict ourselves to events surviving the 
kinematical cut in Eq. o, as in the previous case. Although these two distributions 
also show some sensitivity to the Dirac or Majorana character of the neutrino, it is 
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obvious when compared to Fig. [3(a) that the best results are obtained from the analysis 
of the polar angle of the produced neutrino 



(a) 


(b) 


Figure 8: Dependence of the cross section on the angles tpw (a) and (^e (b), for the SM 
and the SM plus a 1500 GeV Majorana or Dirac neutrino. 


5.4 Measurement of heavy neutrino couplings 


In order to measure the moduli of the heavy neutrino couplings to charged leptons, the 
constants in Eq. (El must be theoretically calculated. This can be done in principle 
with full Monte Carlo simulations taking into account all radiation and hadronisation 
effects, as well as the real detector behaviour. The couplings of the heavy neutrino are 
then 


W2 

^eN 

1/2 


^ ^ ^ 
/le 4^ 4, ’ 



(33) 


The uncertainty in their measurement comes from the statistical fluctuations of the 
signal and background, as well as from the theoretical calculation of the constants Ai. 
Among other factors, their calculation is affected by the inherent uncertainties in the 
signal normalisation. We expect that Ai can be obtained with a precision of 10%, 
which in principle does not affect the determination of coupling ratios. 

We estimate the accuracy with which heavy neutrino couplings could be extracted, 
calculating the constants from the peak cross sections given by our Monte Carlo for a 
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“reference” set of couplings, and assuming a common 10% error in their determination. 
Then, using as input the cross sections for VeN = V^iN = VtN = 0.04, the couplings 
obtained are 


l/ejv = 0.0388 ± 0.00034 (stat) ± 0.0019 (sys), 
Vfj,N/VeN = 1.007 ± 0.016 (stat), 

VrNiVeN = 1-030 ± 0.028 (stat). 


(34) 


In this case the statistical precision of the ratios is very good, of a 2 — 3%. The 
uncertainty in VeN is of a 5%, dominated by systematics. 

A second issue is the determination of the chirality of £NW couplings. This can be 
done with the measurement of the is forward-backward (FB) asymmetry in the W rest 
frame [38,39], which is sensitive to the chiral structure of the iNW coupling involved 
in the decay N i~W^ i~cs (and in its charge conjugate). The measurement of 
this asymmetry requires to distinguish between the two quark jets resulting from 
decay, what can be done restricting ourselves to cs and taking advantage of c 

tagging to require a tagged c jet in the final state. This reduces the cross sections by 
a factor of four. 


We define as the angle between the momenta of the charged lepton i and the s 
jet, in the W rest frame (the definition is the same for N decays into i~W^ —> i~cs 
and i'^W~ —> i'^cs). The FB asymmetry is defined as 


N{cos 9is > 0) — iV(cos 9is < 0) 
N{cos 9£s > 0) -|- N {cos 9es < 0) ’ 


(35) 


with N standing for the number of events. If we parameterise a general iNW vertex 
(ignoring effective terms) as 


CewN = ir {QlPl + 9rPr) NW^ + H.c. (36) 


and assume that the W coupling to quarks is purely left-handed, the FB asymmetry 
is 

_ 3 M ^ - |^^|2 

4M^ + 2m^ + 

In the case of heavy neutrino singlets Ql = VeN, Qr = 0 as seen in Eq. dn, and the 
second factor in Eq. JHH) equals unity. Still, the first factor is very small for ruAr ^ Mw, 
and for mjv = 1500 GeV we have Ape = 4.3 x 10“^ for the three lepton flavours. Such 
asymmetries are unobservable, as long as the experimental statistical errors expected 
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in the experiment are typically AAps ^ 0.012.® Still, for m^r of several hnndreds of 
GeV the asymmetries are measnrable. We analyse snch possibility in next section. 

Further observables are sensitive to the structure of the INW vertices, namely 
spin asymmetries. A proper analysis requires the search for a direction along which 
the N polarisation is maximal, what is beyond the scope of this work. Besides, the 
information on couplings extracted from such observables is expected to be less clean, 
since they involve additional variables apart from masses and couplings. 


6 Heavy non-decoupled neutrinos at ILC 

Heavy neutrinos with masses up to a few hundreds of GeV can also be produced 
at ILG with a GM energy of 500 GeV. Many features of the production process are 
common to both mass and GM energy scales, but in some other respects ILC and 
CLIC are rather different. In order to have a better comparison between both cases 
we summarise here several results for ILC extending the work in Ref. [10]. We take a 
mass rriN = 300 GeV for most of our computations, which follow closely what is done 
for CLIC (see also Ref. [10] for details). The integrated luminosity is assumed to be 
345 fb“^, corresponding to one year of running. 

With the same signal reconstruction method as for CLIC we obtain for m^r = 300 
GeV the combined limits on Ww and or V-iv plotted in Fig. [HI In contrast with the 
behaviour obtained at CLIC, at ILC the sensitivities in the muon and electron channels 
are similar, and both are better than in rWu production. This can be clearly observed 
in both plots: a uNW coupling has little effect on the limits on Ktv, but a coupling 
with the tan decreases the sensitivity, because the decays in the tan channel are harder 
to observe. The direct limit on Wiv, obtained here improves the indirect one only 
for < 0.01. However, as it has already been remarked, the latter is not general 
and can be evaded with cancellations among heavy neutrino contributions [10]. 

The dependence of the total cross section on rriN can be seen in Fig.ITTIRaL 

for VeN = 0.073, = V-w = 0. For a heavier N the cross sections are smaller and 

thus the limits on Wiv are worse. However, up to tuat = 400 GeV this is compensated 

®In the absence of background, the statistical error of an asymmetry A is given by the simple 
expression AA = {I — A^)/{S L), with S the signal cross section and L the integrated luminosity. 
With an asymmetric background the expression is more involved but it can be seen that AA is larger 
than without background. 


23 





(a) 


(b) 


Figure 9: Combined limits obtained at ILC on: and V^at, for Kiv = 0 (a); V^n 

and VrN, for V^at = 0 (b). The red areas represent the 90% CL limits if no signal is 
observed. The white areas extend up to present bounds VeN < 0.073, < 0.098, 

VtN < 0.13, and correspond to the region where a combined statistical significance of 
5a or larger is achieved. The indirect limit from ju — e LFV processes is also shown. 
We take mrq = 300 GeV. 



Figure 10: (a) Cross section for e+e“ —*• at ILC for V^Af = 0.073 and different 

values of mAf. (b) Dependence of the discovery and upper limits on Ww on the heavy 
neutrino mass. Both plots assume mixing only with the electron. 


by the fact that the SM background also decreases for larger mejj. The limits on Ww 
are shown in Fig. ^l(b) as a function of mqq, assuming that the heavy neutrino only 
mixes with the electron. 
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If a heavy neutrino is discovered at ILC, its nature can be uniquely determined 
with the study of the (fjq dependence of the cross section, as it has been discussed 
in section [Ol The dependence of the peak cross section on ipisi is shown in Fig. [TT] 
(a). The results after subtracting the SM contribution can be seen in Fig. [TT](b). The 
distribution is less concentrated at cosipN = ±1 than for CLIC energy but allows to 
determine unambiguously the neutrino character even if a relatively small number of 
signal events is collected. The dependence of the peak cross section on ipw and (pe is 
presented in Fig. 



(a) (b) 


Figure 11: (a) Dependence of the cross section at ILC on the angle (^at, for the SM 
and the SM plus a 300 GeV Majorana (M) or Dirac (D) neutrino, (b) The same, but 
subtracting the SM contribution. 


The measurement at ILC of the V^at, V^tv, K-w moduli can be done in the same 
way described for CLIC, without any substantial differences, thus we do not repeat 
that discussion here. Still, an important difference appears in the determination of 
their chiral structure. In contrast with TeV scale neutrinos, for tuat = 300 GeV the 
FB asymmetry is observable. We show in Fig. [121 (a) the dependence of the SM 
and SM plus heavy neutrino peak cross sections with respect to the angle Oes- The 
measurement of the FB asymmetry requires the subtraction of the SM prediction, 
giving the distributions shown in Fig. [TSl (b). For a Majorana neutrino with Ww = 
0.073, = VrN = 0, the resulting asymmetry is Tps = 0.083±0.016, where the error 

quoted is statistical. The SM and SM±iV cross sections after c tagging are 13.4 and 32.7 
fb, including only W hadronic decays with a c quark. For a Dirac neutrino the results 

®It should be noted that the (pvc and tpe distributions presented in Ref. [10] correspond to the 
whole range of niejj and not only to the peak. 
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(a) (b) 


Figure 12: Dependence of the cross section at ILC on the angles (pw (a) and Pe (b), 
for the SM and the SM plus a 300 GeV Majorana or Dirac neutrino. 


fonnd are eqnivalent within Monte Carlo nncertainty. The FB asymmetry obtained 
from the simulation is in good agreement with the theoretical value Afb = 0.094 
expected for a purely left-handed coupling. 



COS 0 cos 0 

es es 

(a) (b) 

Figure 13: (a) Dependence of the cross section at ILC on the angle 6*es for the SM and 
the SM plus a 300 GeV Majorana neutrino, (b) The same, but subtracting the SM 
contribntion. 

If a heavy neutrino signal is observed in this mass range, we expect that the FB 
asymmetry will have a snfficient statistical signihcance so as to indicate that the eNW 
conpling is left-handed, at least after several years of rnnning. The precision of the 
Mfb measurement strongly depends on the size of the N cross section, therefore the 
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possibility of extracting the left- and right-handed parts of the conpling from the total 
cross section and Apn measurements is difficult to assess in general. 


7 Conclusions 

Neutrinos with masses of the order of 1 TeV are predicted by many models attempting 
to make the new physics scale observable at future colliders. Their presence leads to 
a greater complexity of the models (in order to reproduce the light neutrino masses), 
especially when they are not decoupled from the light fermions. Nevertheless, if they 
exist, large colliders can discover them or provide the best limits on their masses 
and mixings. Such discovery would give a new insight into the mechanism for neutrino 
mass generation. 

In this paper we have concentrated on the production of heavy neutrino singlets 
in association with a light neutrino. If they are the only addition to the SM, the 
production cross section for NN pairs is suppressed by extra mixing angle factors 
This is also reduced by the smaller phase space and by additional decay branching 
ratios, what makes this process much less sensitive to the presence of these heavy 
fermions.^ In left-right models the new gauge interactions allow to produce heavy 
neutrinos in pairs, or in association with charged leptons, with mixing angles of order 
unity. LHC will be sensitive to both mechanisms [40], but the most stringent limits 
are expected from the latter processes [41,42]. On the other hand, colliders will 
be mainly sensitive to neutrino pair production [12,43], which might be used to learn 
about heavy neutrino properties as well [44]. We note that in the case of Majorana 
neutrinos, NN production may give an interesting lepton number violating final state 
signal NN —^ [26,45], which has smaller backgrounds than the analogous 

lepton number conserving final state . In case that new interactions are 

pushed to high energies, heavy neutral leptons can still be produced in pairs if they 
transform non-trivially under the SM gauge group [25]. 

We have discussed the prospects for TeV scale neutrino detection in IWv production 
at future colliders, taking as example the CLIC proposal for a 3 TeV collider at 

^For e+e” ^ NN the t and M-channel exchange diagrams are quadratic in 14v, and the s-channel 
Z diagram has the mixing factor Xatat in Eq. 0 , which is equal to |14Afp + lW^^P + |br^^P- Therefore, 
the cross section is proportional to light-heavy mixing angles to the fourth power. For toat = 1 TeV and 
a coupling VeN = 0.004 in the discovery limit, the extra mixing angle suppression of this cross section 
is of order VJy = 1.6 x 10“®, already giving unobservable rates for heavy neutrino pair production. 
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CERN. This study complements the analysis of the CLIC potential based on new gauge 
interactions [12,43]. We have examined the dependence of the results on the heavy 
neutrino mass, and special emphasis has been made on the importance of a complete 
study of the three possible channels i = e, n, r. As it has been argued, a non-negligible 
coupling to the electron, VeN ~ 10“^ — 10“^ for rriN = 1 — 2 TeV, is necessary to 
produce the heavy neutrino at detectable rates. The produced neutrinos can then 
decay N —> eW and, if they couple to the muon and tau, N —> /ilT and N —>■ tW as 
well (plus additional decays mediated by neutral currents). Among these, the muon 
channel is by far the cleanest one due to its smaller background. We have explicitly 
shown that for uin = 1500 GeV the sensitivity to Ww improves from VeN ~ 0.008 to 
VeN ~ 0.0035 in the presence of a small coupling V^j^n ~ 0.005, but is hardly affected 
by N mixing with the tau lepton. 

We have also studied what we could learn about heavy neutrinos if they were 
discovered at CLIC. It has been shown that the angular distribution of the produced 
neutrino relative to the beam axis is a powerful discriminant between Majorana and 
Dirac neutrinos, giving a clear evidence of their nature even for relatively small signals. 
Then, we have discussed how to extract heavy neutrino charged current couplings from 
data, estimating with an example the precision with which they could be eventually 
measured. We have proposed a method to determine the chiral structure of the iNW 
interactions, which unfortunately is only useful for neutrino masses below the TeV 
scale. 

Finally, we have performed a similar analysis for heavy neutrinos of few hundreds 
of GeV at ILC, extending a recent analysis [10]. In particular, we have explicitly 
investigated the flavour and mass dependence of the limits on neutrino mixing with 
charged leptons. For a heavy neutrino with a mixing large enough, we have shown 
how to establish its Dirac or Majorana nature and how to determine the chirality of 
its charged current couplings. 

The results obtained show that CLIC would outperform previous machines in hnd- 
ing direct or indirect signals of heavy neutrinos. It would extend direct searches up 
to masses around 2.5 TeV, and for masses around 1 TeV it would be sensitive to cou¬ 
plings VeN — 4 X 10“^. If no heavy neutrino signal was found at CLIC, the bound 
VeN < 2 — 6 X 10“^ could be set for m^v = 1 — 2 TeV, improving the present limit 
VeN < 0.073 by more than one order of magnitude and matching a future limit ob¬ 
tained at the GigaZ option of ILC [46], for which a 10^ statistical gain would be 
expected to translate into bounds ~ 30 times more stringent than those of Eqs. la. 
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The direct limit on the product obtained from the non-observation of heavy 

neutrinos would be much more restrictive than the present indirect bound from LFV 
processes, and would remain competitive with future improvements of the upper bounds 
on Br(/r —»• ey) [47] and /i —e conversion in nuclei [48]. For heavy neutrinos with masses 
of few hundreds of GeV, CLIC would probe mixing angles VeN ~ 10“^, one order of 
magnitude better than what will be achieved at ILC. 
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